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I Abstract 

<D We show that the recent Planck Umits on primordial non-Gaussianity impose strong constraints on 

i__i light hidden sector fields coupled to the inflaton via operators suppressed by a high mass scale A. 

I We study a simple effective field theory in which a hidden sector field is coupled to a shift-symmetric 

K>" inflaton via arbitrary operators up to dimension flve. Self-interactions in the hidden sector lead 

P^ to non-Gaussianity in the curvature perturbations. To be consistent with the Planck limit on local 

CN non-Gaussianity, the coupling to any hidden sector with light flelds and natural cubic couplings must 
be suppressed by a very high scale A > 10 // . Even if the hidden sector has Gaussian correlations, 

^^ nonlinearities in the mixing with the inflaton still lead to non-Gaussian curvature perturbations. 

^^ In this case, the non-Gaussianity is of the equilateral or orthogonal type, and the Planck data 

• • requires A > 10 H. 
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1 Introduction 

The recent results of the Planck satellite [1-3] represent a new standard in precision cosmology. 
With cosmic variance limited measurements of the temperature anisotropy up to £max ~ 1500, the 
theory of initial conditions is tested at a precision of -^^ax ~ 7 x 10~^. This has led to a significant 
improvement in the constraints on the non-Gaussianity of the primordial fluctuations [3]. These are 
reported as constraints on the amplitudes of three bispectrum templates — local [4] , equilateral [5] , 
and orthogonal [6]: 

/ir' = 2.7±5.8, /T'' = -42±75 and /°[*^° = -25 ± 39 . (1.1) 

Recalling that a rough quantitative measure of the amount of non-Gaussianity is |/nl|A(^ < 10^^, 
the limits in (1.1) represent a remarkable level of Gaussianity in the primordial perturbations. 

It is very well known that limits on non-Gaussianity lead to constraints on the self-interactions of 
the inflaton.^ However, the requirement that the inflaton generates prolonged accelerated expansion 
entails further, model-dependent constraints on inflaton self-interactions, which can be stronger 
than the constraint from (non) observation of the bispectrum. For example, in slow-roll models the 
inflaton $ respects an approximate shift symmetry that constrains self-interactions in the potential 
to be small and non-Gaussianity to be unobservable [9]. In contrast, self-interactions of additional 
fields S that couple weakly to the inflaton are much less constrained: S could have a large cubic 
coupling^ 

C^ 3 -/uS3 , (1.2) 

with /i ~ //, and hence have non-Gaussian correlations, without interfering with the slow-roll evo- 
lution. 
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inflaton sector '^ tiidden sector 

Figure 1: The inflaton sector and the hidden sector mix through irrelevant operators suppressed by 
the scale A. 

In this paper, we ask how the Gaussianity of the CMB can be used to place constraints on hidden 
sector fields. We consider a hidden sector field S that is coupled to a shift-symmetric inflaton <i> via 
irrelevant operators suppressed by a high scale A. We construct a general effective field theory (EFT) 
involving these fields, 

/:eff [$, s] = a + £s + /:mix[^, n ■ (1.3) 

The leading mixing respecting the shift symmetry is the dimension-five operator 



^In terms of the cutoff scale A of an effective field theory for the fluctuations of a single inflaton [7], the constraints 
on /nl"' and /5l''° in (1-1) imply that A > 'iH [8], where H is the inflationary expansion rate. 

^The conformal symmetry of the inflationary quasi-de Sitter background fixes the three-point function of E up to 
an overall normalization [10]. We therefore lose no generality by focusing on E'^. 



Such a coupling to the inflaton kinetic term arises rather naturally in ultraviolet (UV) completions 
of inflation [11]. If S is light enough'^ to be quantum- mechanically active during inflation (m < ^H), 
self-interactions in the S sector can be imprinted as non-Gaussianities in the visible curvature per- 
turbations, through the mixing (1.4). In fact, the theory has two distinct sources of non-Gaussianity: 
self-interactions in the hidden sector, (1.2), and nonlinear couplings between the two sectors, (1.4). 
As a result, the phenomenology of the model is rather rich, with all three bispectrum shapes probed 
by Planck realized in different regions of the parameter space. We will use the Planck bounds on 
non-Gaussianity to put a lower bound on the scale A in (1.4), 



A > O{W^-10^)H , (1.5) 



where the precise numerical coefficient on the r.h.s. depends on the couplings of the hidden sector. 
A detection of primordial tensors would show that H > lO^^Mpi, so that A can exceed the Planck 
scale. In this case, the limits on the bispectrum [3] constrain the ultraviolet completion of gravity. 

A similar strategy was first outlined in [12], for particular classes of strongly-interacting hidden 
sectors. In the analysis of [12], the bispectrum could be computed only when the mixing between the 
sectors was perturbatively small, so that a region of parameter space was excluded not by data but 
by the need for theoretical control. In the present setting, we will be able to compute the observable 
signatures for all values of A. 

Our approach is analogous to the use of high-precision measurements in particle physics to con- 
strain the scale A of physics beyond the Standard Model (SM). There one can place lower bounds on 
A that depend sensitively on the effective operators under consideration, ranging from A > 5 x 10^ 
GeV for precision electroweak measurements [13] to A > 10^^ GeV for proton decay [14]. Moreover, 
the discovery that neutrinos have finite masses hints at new physics at A ~ 10^^ GeV.*^ In this paper, 
we apply the same philosophy to inflation and use the Planck data to put precision constraints on 
the scale A of couplings to hidden sectors. 

Outline. — In Section 2, we construct an effective field theory that couples a shift-symmetric inflaton 
fleld to a hidden sector field with strong self- interactions. We derive the effective action for small 
fluctuations around a slow-roll background, and we briefly discuss the radiative stability of different 
parameter regimes. In Section 3, we use both analytical and numerical techniques to study the 
phenomenology of this theory. By scanning the complete parameter space of the EFT, we discover 
that all three bispectrum shapes probed by Planck are realized. We then show how the Planck 
limits on local, equilateral and orthogonal non-Gaussianity probe physics at energy scales orders of 
magnitude above the inflationary expansion rate. In Section 4, we consider the illustrative example 
of a supersymmetric hidden sector. We explain how the hidden sector self-interactions, and the 
corresponding non-Gaussianity, depend on how supersymmetry breaking is communicated to S. 
Our conclusions appear in Section 5. Appendix A contains details of the numerical analysis, while 
Appendix B provides further discussion of the EFT of the two-scalar system. 

Throughout, we will use natural units with h = c= 1 and reduced Planck mass Mpi = l/\/8nG. 
Our metric signature is ( — \- ++) and we deflne (d^)'^ = g'^'^df^^d^^. 



^The presence of hidden sector fields with masses m < H is not unreasonable from the perspective of UV completions 
of inflation [11]: for example, E could be a modulus that is massless before supersymmetry breaking. 

*The EFT of the SM contains a unique dimension-five operator involving two leptons and the Higgs, £|,g — 
{LH){LH)/A. For A ~ 10^^ GeV, this operator induces neutrino masses of interesting size, m^ ~ 10~^ eV. 



2 Effective Theory of Multi-Field Inflation 

2.1 EFT of Background Fields 

Our interest in this paper is the coupUng of an approximately shift-symmetric inflaton field <I> to a 
hidden-sector field S, which may have significant self-interactions. Mixing operators suppressed by a 
high scale A turn the non-Gaussianity in the hidden sector into observable curvature perturbations. 
In this section, we discuss the effective field theory of $-S couplings.^ 

We assume that the quantum-corrected inflaton action takes the slow-roll form 

U = -^-{d^f-V{^). (2.1) 

The observed (near) scale-invariance of the primordial fluctuations [1, 2] suggests that $ respects 
an approximate shift symmetry, <1> i— )• $ + const. This constrains self-interactions in the potential 
to be small and non-Gaussianity to be unobservable [9]. Large interactions could still come from 
higher-derivative operators such as {d^)^/A^ [19]. However, these operators only produce observable 
non-Gaussianity if A < ($)^'^, and in that case the derivative expansion of the EFT cannot be 
truncated, so that the theory requires a UV completion to be predictive [20].® Although it would be 
straightforward to generalize our treatment to include higher-derivative kinetic terms, in this paper 
we will consider only non-Gaussianity from single-derivative couplings. 
We take the hidden sector to have the Lagrangian 

/:s = -2(9S)2-y(S). (2.2) 

Interactions in the hidden sector are much less constrained than in the inflaton sector, and in 
particular the self-interactions in the potential V^(S) can be large. At cubic order in S, the self- 
interactions of S are completely characterized by /iS^: higher-derivative operators can be written 
in terms of S'^ and total derivatives (see footnote 2). When /x ~ i?^, the fluctuations of T, have 
order-one non-Gaussianity. 

We characterize the communication between the two sectors in terms of the most general mixing 
operators allowed in the EFT, 

C^,^[<^,^ = J2ci^j^, (2.3) 

where the operators Oj are made out of powers of the fields <^ and S and their derivatives. The 
parameters 5j are the mass dimensions of the operators O/, and the c/ are dimensionless Wilson 
coefficients. Before imposing that $ enjoys an approximate shift symmetry, the leading terms in 
£niix are simply the lowest-dimension operators involving both $ and S. In Table 1 we list all 
possible operators constructed from at most first derivatives, up to dimension 6 = 5. 

The requirement of an approximate shift symmetry' $ i— ;■ <I> + const, excludes all operators 
involving $ (as opposed to 9/i<&). Three mixing operators remain: 



^See [15-17] for an alternative perspective using the effective theory of inflationary fluctuations [7]. We note that 
the derivative expansion in the effective theory of the background fields is different from that of the fluctuations [18]. 

®In this case, the constraints on non-Gaussianity are more usefully organized in terms of an EFT for inflationary 
fluctuations [7]. The advantage of this approach is that it captures all single-field models, including those with 
A < (<i?)^''^, in the regime of control of a single EFT. 

'^We do not have to commit to whether this symmetry is a fundamental symmetry of the UV theory, an accidental 
symmetry of the IR description, or simply a consequence of fine-tuning in the quantum corrected effective potential. 



Dimension Operators 



Vq 



$ , s 



$2 , I]2 ^ $E 



$3 , 1]3 ^ $25. ^5^2 



$4 , I]4 $3j^ ^25.2 ^j. 



a,,$(9^$ , a,,s9^E , 5,, $9'^!] 




Table 1: List of all operators constructed from at most first derivatives, up to dimension 5. 

This dimension-four kinetic mixing can be removed by a rigid rotation in field space, <i> 1— )• $ = 
cos(6i)$ - sin((9)S and S h^ S = sin(6l)$ + cos(6l)S. But as a result ^(S) = V{t, ^), so that 
<& is not shift-symmetric. In order for the canonical field $ (with diagonal kinetic term) to 
be shift symmetric, we will assume that the coefficient of (9^$5^E is negligibly small. This is 
equivalent to imposing the shift symmetry after diagonalizing the kinetic term. We emphasize 
that this shift symmetry is motivated by the fact that the primordial fluctuations are observed 
to be nearly scale-invariant. 

• 6 = 5: {d^.^df'^)^ 

The dimension-five coupling (5^$9'^S)S can be integrated by parts to give^ 

Ct.^{d,^d^^)j ^ -l^^'- (2-4) 

This operator is redundant: using the $ equation of motion, □<!> = d^V, we can rewrite it in 
terms of ^"'S"'. 

• 6 = 5 : {df,^d''^)E 

A primary goal of this paper is to characterize the effects of the interaction 

**For equal-time correlation functions, it may not be obvious that a total derivative can be neglected, due to the 
boundary term. In the in-in formalism, these boundary terms are always associated with equal-time commutators and 
(assuming locality) can be removed by a redefinition of the local operators. 



This operator is the dominant source of mixing between fluctuations in the hidden sector and 
the visible sector. 

The operator (2.5) contains a tadpole, 

that drives the field S away from the origin. We will assume that the potential VC^) in (2.2) 
stabilizes S at Sq = const, (in Appendix B, we examine the naturalness of this assumption). Once 
S has developed a vev, the interaction (2.5) induces a correction to the inflaton kinetic term, 

U = -l{l + K){d^f-V{^), (2.7) 

where k = Sq/A. We will usually work in the regime where this is a small effect (i.e. k <^ 1). Even 
when K ~ 1, this correction can be absorbed by a field redefinition, ^^ = (1 + k)<1>^. In terms of the 
new field, the mixing term becomes 

-^mix — "2 A ' '^ -^ 

where A = (1 + k)A. From now on we will assume that the rescaling has been performed and drop 
the tildes. Moreover, we assume that the potential for $ (really $) satisfies the slow-roll conditions, 
and that the inflaton has a time-dependent background vev ^o{t). 

2.2 EFT of Coupled Fluctuations 

Cosmological observables are sourced by the small fluctuations around the background vevs, i.e. 

^{t,x) = ^o{t) + ^{t,x) and T.{t,x) = T,o + a{t,x) . (2.9) 

We will work in spatially flat gauge, gij = a^dij, so that the primordial curvature perturbation is [21] 

C{t,x) = -^^{t,x) . (2.10) 

At leading order in the slow-roll expansion, the inflaton fluctuations if are massless^ and the mixing 
between matter and metric fluctuations vanishes. Eqs. (2.1) and (2.5) then become^*^ 

U + CBCgi, = -lid^f + P'pa - ^^^^ , (2.11) 

where we have deflned the important mixing parameter 

p^^. (2.12) 



^Expressed in terms of (" this statement is exact. 

^"We have dropped a cosmological constant term and two tadpoles proportional to (p and a. These terms are 
cancelled by potential terms in eqs. (2.1) and (2.2). This is discussed further in Appendix B. 



We have absorbed the Wilson coefficient cb into the definition of A. Similarly, eqs. (2.2) and (2.4) 
give^""^ 

C^ + CACt. = -lidaf - ^m^a^ - ^la^ + ■ ■ ■ , (2.13) 

where m^ = Vq — 3ca (1 — jt^P^ ^^"^ M — W.^o" '^'^^^ (dropping contributions that are suppressed 
by slow- roll parameters). Up to cubic order, the complete Lagrangian for the coupled ip-a system is 
then 

CM^,a] = -\{d^? -\{daf -\m\' + p^a - ^^^"/^^^ • (2-14) 

^ V ' ^ V ' 

•Co Ant 

In the following, we will study the phenomenology of the Lagrangian (2.14), both analytically and 
numerically. 

2.3 Dynamics and Naturalness 

The effective Lagrangian (2.14) contains three independent parameters: the mass m and the cubic 
coupling /i of the hidden sector field, and the strength p of the mixing with the visible sector. 
Depending on the size of the mixing parameter p relative to H, we encounter different dynamical 
regimes. Throughout our analysis, it will be useful to consider strong mixing {p ^ H) and weak 
mixing {p < H) separately. These two cases are qualitatively different, in terms of dynamics, 
observational signatures and the natural range of parameters. We can understand the origin of this 
difference directly from the quadratic Lagrangian, 

^2 = -lidiff - ^{daf - \m^a^ + pupa . (2.15) 

The equations of motion associated with this Lagrangian are 

/j2 

(p + 3H^ + ^^ = -p\& + 3Ha] , (2.16) 

a + 3Ha+(^+mAa = pip, (2.17) 

where we have assumed a flat FRW background with scale factor a{t) and Hubble parameter H = 
a I a. These equations of motion admit WKB-like solutions, where (/? oc o" oc e~*J '^^* ' , as long as 
uj^H. By comparing terms in (2.15), we see that when uj ^ p the mixing term, /O^ia, is negligible 
compared to ip^ and d^ and may therefore be treated as a perturbation. On the other hand, when 
uj <^ p, the converse is true and we may neglect the kinetic terms 99^ and a'^ and describe the 
dynamics in terms of the mixing term and the gradients. Which of these terms dominates at horizon 
crossing, a; ~ i/, is a key distinguishing feature of the weak and strong mixing cases, respectively. 



^The operator £„ 



A 1 n^o 2 1 □</:'. 2 . „^ N 



Using the equations of motion, n$o ~ SH^o = ~3pHA and Oip — p{a + 3Ha), these terms can be absorbed into 
the definitions of m^ and p. 



2.3.1 Weak Mixing 

In the case of weak mixing, p < H, the dominant kinetic terms are the usual ones, ip^ and a"^. As 
such, the dynamics follows a familiar pattern: the relevant scales are the mass m and cubic coupling 
//, as well as the mixing scale p. The hierarchy of these scales is illustrated in figure 2. In order 
for the fluctuations of a to be generated during inflation (and to be weakly coupled) we will require 
that {m, p} ^ H. As a result, the full parameter space of this regime is {m, p, p} < H. 

CO 
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m 
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Figure 2: Hierarchies of scales in the weak mixing regime. 



Let us briefly discuss whether these conditions are stable under quantum corrections [15]. First, we 
estimate the loop corrections to the mass parameter m?. The cubic interaction is a soft deformation 
of the theory and corrects the mass by a finite amount. 



6m ~ p 



(2.18) 



For p < H, this is an acceptably large correction to the mass. On the other hand, the loop correction 
associated with the operator j^{dip)'^a is 



dm 



A2 



(2.19) 



where A^. is the UV cutoff of the effective theory of the fiuctuations. We see that the natural value of 
m^ is sensitive to A^,. A natural choice for A^ is the symmetry breaking scale (^o)^ [7, 8], i.e. the 
scale at which we integrate out the background vevs and focus only on the fiuctuations. Letting 
AJ ~ <i>o, we get 



A2 



Similarly, we can estimate the one-loop correction to the cubic coupling 

5p 



A| '^2 

A3 A3 



'1 

A 



^^^C«H. 



(2.20) 



(2.21) 



We see that, as long as the cutoff of the EFT of the fluctuations is not far above the symmetry 
breaking scale, loop corrections do not destabilize the parameters of the EFT. 



2.3.2 Strong Mixing 

In the regime of strong mixing, p ^ H, the dynamics is controlled by a linear kinetic term pipa. 
When this kinetic term dominates, the scaling behavior of operators is modified and the kinetic 
terms ip'^ and o"^ become irrelevant. The quadratic action can be approximated as 



£5 



pipa 



2 a2 



1(^ 
2 a2 



1 



2 2 
m a 



(2.22) 



The WKB-like solutions behave as a; ~ k'^/p at early times. As a result, if we assign the time 
coordinate, t, a scaling dimension +1, then the space coordinate, x, will have scaling dimension j- 



In order to make this scaling manifest, we define x* 
of these new variables, the action becomes 



,1/2 



X*, (f 



p ^'^ip and a = p ^'^a. In terms 



S 



dtd-^x a^ 



ipa 



1 my 

2 a2 



2 a2 



-Ma^ 



p^l^d^ 



~r.\-2 



1 1 (dj^) 

2 A3/4 ■ 



(2.23) 



where M = m?/p, p 



1/4 



MP' 



'^'^ and A^'^ = Ap ^'^. We have dropped the (p'^a interaction because 



it is suppressed by additional powers of p. In writing this action, we have assumed that the mass 
term is negligible at the scale p. This assumption requires that M <^ p. In addition, for the theory 
to be perturbative at horizon crossing, we also need p ^ H and A^ H. 
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Figure 3: Hierarchies of scales in tlie strong mixing regime. The Hubble scale H may be above or 
below the scale M = m^ /p. 



Despite its appearance, this Lagrangian describes a single degree of freedom. The field a now 
plays the role of the conjugate momentum of if, i.e. ptp = pa. When u < M, we can integrate out a 
(e.g. by completing the square [8]) to produce an effective Lagrangian for ip with dispersion relation 
uj = Cgk, where Cg = M/ p. For this reason, the phenomenology of this model will depend sensitively 
on M/H, as illustrated in fig. 3. 

Now that we understand the allowed range of parameters, we should discuss the natural values 
of these parameters. The effective theory described by (2.23) is well-defined up to the energy scale 
p where we are required to include the second degree of freedom. We assess naturalness of this EFT 
by computing loops with the UV cutoff A^, < p. At one loop, there are three contributions to the 



mass parameter M: 

-1/4 A 3/2 /. X 3/2 

5'^^^M ~ fi^'^Kl'^ , 5(^'M ~ ^ , '^ , ,5(^'M ~ Aj ^ J . (2.24) 

In writing these expressions, we have neglected additional numerical suppressions such as (IGvr^)^^. 
Regardless, we see that jl < H and A ^ p is sufficient to ensure that M <^ p. On the other 
hand, whether M is above or below the Hubble scale depends sensitively on details, including these 
numerical factors. 

The parameters in the effective Lagrangian are also renormalized by loops above the scale p, 
including the contributions computed in the previous section. One might worry that for A^ ~ <i>o, 
eq. (2.19) implies M ~ p. These loop corrections could be controlled if the theory is supersymmetric 
above the scale p, which would have no impact on the low energy EFT [15]. Alternatively, including 
the appropriate factor of (IGvr^)^^ can be sufficient to create a hierarchy between p and M. 

3 Non-Gaussian Phenomenology 

The Lagrangian (2.14) contains two distinct sources of non-Gaussianity: 

4;; = 4^ -d £» = -..'. (3.1) 



These interactions get converted into a bispectrum of primordial curvature perturbations via the 
quadratic mixing term £mix = p^cr. In this section, we will compute the bispectrum as a function 
of the parameters m, p and p (or A). Our treatment will be numerical, except in the limits of weak 
mixing (p < H) and strong mixing (p > H), where we present analytical results. The constraints on 
non-Gaussianity from the Planck satellite can be viewed either as upper limits on p or lower limits 
on A. The second point of view makes the Planck experiment a probe of high-scale physics. 

3.1 Preliminaries 

For the convenience of the reader, we collect a few basic formulae that are used in the statistical 
analysis of non- Gaussian perturbations. 

The main diagnostic for primordial non-Gaussianity is the bispectrum, 

iCk.Ck.Cks) = (2^)' B^{kuk2, h) <5(fci + k2 + fcs) • (3.2) 

We will use the in-in formalism [22] to compute the bispectrum for the Lagrangian (2.14), 



(CfclCfc2Cfc3)(*) 



Te 



i Jl dt' Hintit') 



CkiCk2Cka{t) 



Te 



i Jl dt' Hintit') 



|0) , (3.3) 



where |0) is the vacuum of the free theory^^ and .ffint[C5<5'] is the interaction Hamiltonian. The 
operators (^ and a are expanded in terms of creation and annihilation operators, 



^^The initial condition will include a small evolution into the imaginary time direction that projects the theory onto 
the interacting vacuum. 



10 



and similarly for 6"^ . The mode functions Cfc {t) oc Vfc (^) ^^'^ ^k (i) are interaction picture fields whose 
time evolution is determined by the quadratic part of the Hamiltonian, cf. eqs. (2.16) and (2.17). 
The bispectrum is computed perturbatively as an expansion in Hint- Time integrals are performed 
by analytical continuation into the complex plane (see Appendix A). We will compare our results to 
the three template bispectra used in the Planck analysis [3], 

Aocai = ^ [PiP2 + 2 perms.) , (3.5) 

^cqui = 1{q {PfPiPz)^'^ - 3P1P2 - 2 (P1P2P3)'/' + 5 perms.) , (3.6) 

^ortho = I (18 {PfPiP^f^ - 9P1P2 - 8 (P1P2P3)'/' + 5 perms.) , (3.7) 

where Pi = Pc^{ki) is the power spectrum, (CkiCkj) = (2vr)'^ Pc,{ki) 6{ki+kj). To quantify the degree of 
correlation between two bispectrum shapes it is convenient to define the following shape function [5] : 

S{xi,X2) = ixiX2f Bi;{xi,X2,l) , Xi = ki/k3, (3.8) 

with inner product 

F{S,S')= / S{xi,X2)S'{xi,X2)dxidx2 , (3.9) 

Jv 

where the integrals are only over physical momenta satisfying < xi < 1 and 1 — xi < X2 < 1. Two 

shapes are highly correlated if their normalized scalar product or 'cosine' is close to unity, 

C{S, S') = F{S,S') ^ _ 

y/F{S,S)F{S',S') 

To facilitate comparison with the constraints in eq. (1.1), we also compute the amplitude of the 

bispectrum, 

_ 5 B^{k,k,k) 

^"""18 P^{k) • ^"^-^^^ 

3.2 Spectrum of Non-Gaussianities 

In Appendix A, we describe in detail our approach to computing the bispectrum numerically. Both 
the UV and the IR have to be treated carefully. In the UV, it is important to make sure that the 
quantization of the coupled fields is consistent with the equations of motion (2.16) and (2.17). In 
the IR, spurious divergences can appear in the numerical evaluation of the integrals in (3.3). The 
interested reader is referred to the appendix for a detailed description of how we deal with these 
technical issues. 

Figure 4 shows the result of our analysis for the amplitude of the bispectrum, /nl. For the purpose 
of illustration, the mass of S has been fixed to m = H. The plot then shows the dependence of 
/nl on the parameters p and fi. The grey regions denote |/nl| > 10 and are therefore disfavored by 
the Planck data.^'^ As we will show, the asymptotic limits of this plot (p ^ H and p <^ H) can be 
understood analytically. 



'To be precise about this would require a dedicated likelihood analysis, which is beyond the scope of this paper. 
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/nl = -10 



^i/H 



-10 



Figure 4: Non-Gaussian phenomenology, for m = H. The grey regions correspond to |/nl1 > 10. 



3.2.1 Weak Mixing 

In the regime of weak mixing, p < H (or A > ^q/H), the phenomenology of the (p-a-system was 
studied in detail in the context of quasi-single-field inflation [23] (see also [15]). We will therefore 
cite the main results without proof. 

Mode functions. — Treating the coupling pipa as a perturbation, the mode functions in the interaction 
picture are determined by the free field Lagrangian Cq, 

if" + 2nip' + A;V = , (3.12) 

(t" + 2'Ho-'+ (fc2 + aV2)cj = . (3.13) 

In a quasi-de Sitter background, a^r) w —1/{Ht), the Bunch-Davies solutions to (3.12) and (3.13) 
are 



o-fc(r) 



H 



2A;3 
^ H 



(1 + ikT)e-'''^ , 

(-A;r)3/2HW(-A;r), 



2V2F 



4~lP 



(3.14) 
(3.15) 
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where r is conformal time and Hu is the Hankel function of the first kind. While ip freezes for 
|fer| <C 1 (up to slow-roll corrections), a decays with a rate determined by m/H. The mixing pipa 
is treated as part of the interaction Hamiltonian and perturbatively converts the fluctuations a into 
If (and hence (). 



if. 



a 

-x X- 



f 




Figure 5: Feynman diagrams describing the conversion of hidden sector fluctuations a into observable 
fluctuations ip (or ("). 

Power spectrum. — The leading-order contribution to the power spectrum comes from the inflaton 
fluctuations, while the mixing with a appears in subleading corrections (see fig. 5): 



27rA/ 



6o 



1 + c{m) 



m 



(3.16) 



where an explicit expression for the function c{m) can be found in [23], but will not be needed in 
the following. For the calculation to be under perturbative control we require c{m)p^ < H^. Unless 
m -^ H, this is equivalent to the weak mixing condition. For m <ti H, the function c{m) becomes 
large and the condition for perturbative control becomes the stronger one. 

Bispectrum. — Since the inflaton fluctuations are Gaussian, the leading contribution to the bispectrum 
comes from the hidden sector interactions. Each leg in the interaction £|^[ = — /io"^ is converted to 
C via an insertion of pipa. Using eq. (3.3) to evaluate the cubic diagram in fig. 5, one finds 



/nl = fim) 



1 



^ /P 



27rA<; H \H 



(3.17) 



where the function f{m) can be found in [23]. Eq. (3.17) explains the \p\ < H part of fig. 4. 

Depending on the mass of the extra field, the shape of the non-Gaussianity interpolates between 
local (m < H) and equilateral {m ~ H). This feature is easy to understand intuitively. When 
m ^ H, the field a decays rapidly on superhorizon scales. The interactions between different modes 
are therefore suppressed unless they exit the horizon at nearly the same time. This leads to the 
equilateral shape. On the other hand, if a is nearly massless then modes that exit the horizon 
at different times in the infiationary history still have non-trivial overlap, and the shape of the 
bispectrum will be approximately local. 

Another important aspect of eq. (3.17) is the large prefactor A7 , which allows for observable 
levels of non-Gaussianity without violating weak coupling, i.e for p < H and p < H. The A7 
enhancement also admits a straightforward explanation. For p ^ H, order-one non-Gaussianity 
in the cr-sector (i.e. p/H ~ 1) should lead to order-one non-Gaussianity in the C-sector. But the 
relevant measure of the non-Gaussianity of the bispectrum of C, is /^lA^. As a result, we should 
have /nlA^- ~ {p/H){p/H)'^, which is confirmed by the direct computation of (3.17). 
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3.2.2 Strong Mixing 

When jj -^ H, the dominant source of non-Gaussianity switches from £|^^ to C\^^^ = —i^{dip)^a. 
Observable non-Gaussianity is now generated only if A is small enough so that A < ^q/H (and 
hence p > H), which is the regime of strong mixing. The dynamics in this region of parameter space 
is less well understood, although some aspects have been studied in [8] (see also [24-31]). Here, 
we provide a complete understanding of the phenomenology in the strong mixing limit, as well as 
numerical results for arbitrary values of the parameters. 

Mode functions. — As we have discussed in §2.3, for energies below the mixing scale p, the system 
reduces to a single effective degree of freedom. From the Lagrangian (2.22), we find the equation of 
motion for if, 



where we have defined Ti = a' /a, as well as 



^' + a'^l{l + ^]^ = 0, (3.18) 



n 



^k{r) = ~ and M=— . (3.19) 

p a^ p 



The solution for a is determined by the constraint equation 

aa = 7^^^ . (3.20) 

We see that a is not an independent degree of freedom. In §2.3, we also explained that the parameter 
M = rr? I p is the more relevant mass parameter in the strong mixing regime. For general M, 
eq. (3.18) can only be solved numerically (see Appendix A). However, in the limits M ^ H and 
M <C -f^ an analytical understanding is possible: 

• For M ^ H, eq. (3.18) becomes 

If" + 2nip' + — A;V = , (3.21) 

P 

where we have used that near horizon crossing $7^ ~ i? ^ M. For p > M, we can identify 
(3.21) as the equation of motion of a scalar with non-trivial sound speed, 

c? = ^. (3.22) 

P 

The mode function corresponding to the Bunch-Davies initial state is 

Mr) = ^^(1 + ^c.fcr)e--»^- . (3.23) 

It is easy to show that this identification with a small-Cs theory extends to the cubic Lagrangian 
(see [8, 24]). Both the power spectrum and the bispectrum will therefore be those of a single- 
field theory with reduced sound speed, Cg < 1 [32]. 
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For M <^ H, eq. (3.18) becomes 



^4 

i^" + 4n^' + -^^^ = . 



In de Sitter space, this has the exact solution 






fc^^o 2r(|) H fp\^/^ 



vr 



A:3/2 



(I) 



(3.24) 



(3.25) 



We see that the inflaton fluctuations are enhanced by a factor of {p/HY'^ relative to the 
canonical slow-roll result. 

Power spectrum. — Given a solution for the mode functions, it is easy to compute the power spectrum 
of curvature perturbations in the super horizon limit. 



2tt/\c = d{M) X 



H"^ ( P\ 1/4 



$0 ^H 



(3.26) 



where the function d{M) is computed numerically and plotted in fig. 6. We have confirmed that our 
numerical result has the expected asymptotic limit d{M) — )■ {H/M)^''^ for M ^ H}^ 
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Figure 6: Numerical computation of d{M) (dotted), fi{M) (dashed) and fiiM) (solid). Notice that 
/2(M) changes sign near M = 0.8H. 



Bispectrum from {dip)'^a. — First, we consider the bispectrum generated by the interaction C 



(1) 



hr{dip)'^a. The shape is equilateral and the amplitude is 



2A 



requil 



fZa^ = MM)x 



H 



(3.27) 



The function fi{M) is computed numerically and plotted in fig. 6. Recall that M = m?/p, so for 
fixed m, the function /i contains a dependence on the mixing parameter p. We have confirmed that 
our numerical result has the expected asymptotic limit /i(M) — t- —^H/M for M ^ H}"^ In (3.27), 



-1/2 



iC oc Cs = (p/M) 

rcquil 1-2 _ 1 



1/4 



^*This follows from A;; ex c. 

15 This follows from /°?,"" ^ Icj'^ = \p/M 
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we have factored out p/H to match the scahng in the M 
from (2.23) and dimensional analysis, 



limit. We can understand this scaling 






equil .„ . i/-??\3/4 p 

NL,(i) ~ (2vrA^) [^) = -^ , 



(3.28) 



where we have defined A^' ^ = Kp ^' ^ as before. 

Bispectrum from a^ . — A similar analysis for the bispectrum coming from C-^^^ 

-3/4 



(2) 



-pa gives 



j^equil 

/nL,(2) 



/2(M) p 

X 



27rA/ 



H \H 



P 



(3.29) 



where /2(M) is plotted in fig. 6. The numerical result has the correct asymptotic limit /2(M) — )• 
(H/M)^''^ for M ^ H. Again, we have factored out the scaling behavior in the M — )• limit, which 
again follows from (2.23), 



/-■i~(2.A,)-'(|)"' = (2.A,r""''^-'" 



(3.30) 



where p^''^ = pp~^'^. 

Interestingly, /2(M) — and hence /nl"(2) — changes sign near M/H = 0.8. We interpret this as 
a cancellation between two large equilateral contributions. But such a cancellation was precisely 
what gave rise to the orthogonal shape in [6], which motivates looking more closely at the shape of 
the bispectrum near M/H = 0.8. We have computed the cosine (3.10) between the shape arising 
from the operator Cf^^ = —pa^ and both the equilateral and the orthogonal templates. The result is 
shown in fig. 7. We see that the theory indeed realizes both the equilateral and the orthogonal shapes 
as a function of M = rr? j p. Since the theory also includes the possibility of local non-Gaussianity 
in the weak mixing limit, our simple model realizes all three standard templates probed by Planck. 
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Figure 7: The cosine of tlie sliape S(2) witli the orthogonal template shape So 



Complete bispectrum,. — To summarize, the complete bispectrum predicted by the Lagrangian (2.14) 
in the limit of strong mixing is 



/T" = /i(M).^ + /2(M) 



P_ 
H 



1 



P 



27rA^ H \H 



-3/4 



(3.31) 
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We note that p and // can have either sign, so neither the sign of the total /nl, nor the signs of the 
individual contributions are fixed. Eq. (3.31) explains the \p\ > H part of fig. 4. 

The result may be written more symmetrically in terms of the parameters in (2.23), 



ii'\3/4 



(27rAc)/r" = /i(Af)xd(M).(^) 



+ /2(M)(| 



1/4 



(3.32) 



This makes it clear that the scaling with energy H is given by 4 — Aj, where Aj is the scaling 
dimension of the operator that produces the three-point function (see §2.3.2). 

When p, and p have opposite signs, the two contributions in (3.31) can cancel against each other. 
Naively, this might suggest that we can have large p. and p without producing large non-Gaussianity 
(see fig. 4). This is a bit misleading: one should recall that the parameter /nl is defined in the 
equilateral configuration ki = k2 = k^. The two interactions that we are cancelling against each 
other produce similar, but not identical, shapes. Away from the equilateral limit, the cancellation will 
therefore not be perfect. This leads to a non-trivial bispectrum that is constrained by observations, 
although the /nl parameter suggests otherwise. Again, we can diagnose this by computing the cosine 
with the orthogonal template. The result is shown in fig. 8: we find an orthogonal component near 
the point in parameter space where the two equilateral shapes cancel. 



(l) + (2) 




Figure 8: The cosine of 5(1) + 5'(2) with the orthogonal template shape Sortho, for p — 50H. 



3.3 Probing High-Scale Physics 

Fig. 9 shows the non-Gaussian predictions of the effective theory expressed in terms of the suppression 
scale A of the dimension-five operator (2.5), cf. (2.12). Again we show the region |/nl| > 10 for the 
case m = H. This indicates how the Planck limits (1.1) constrain high-scale physics. 

Non-Gaussian hidden sector. — We get the strongest constraint if the hidden sector is order-one non- 
Gaussian {\p\ ~ H). In the weak mixing limit, we use (3.16) to replace $o in (2.12) and write the 
bispectrum amplitude (3.17) as 



/n 



/(m) I p ( H 



27r A^ iJ V 27rA 



(3.33) 
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Figure 9: Constraints from Planclt on tiie scale of new piiysics. 



Solving for A, we obtain 



A 



9x 10^ 



I/n 



10 



-1/3 



H 



1/3 



H . 



(3.34) 



where the precise numerical coefficient depends on m. For m < H, the scale A is bounded by 
Planck's limit on local non-Gaussianity [3], l/^if | < 10, while for m ~ iJ we should impose the 
limit on equilateral non-Gaussianity, 1/^2"' | < 75. A conservative bound on the scale of new physics 
is therefore 

A > 0(5) X 10^ (M^ H . (3.35) 

For fj, r^ H, this constrains physics at scales many orders of magnitude above the Hubble scale. 

Planck-suppressed operators. — To express (3.35) as a limit on A in terms of an absolute energy scale, 
we need additional observational input. Specifically, a detection of primordial tensor fluctuations 
would relate the inflationary expansion rate H to the Planck scale Mpi. In terms of the tensor-to- 
scalar ratio r, we can write 



and (3.35) becomes 



H = TTAr 



A > 0.5 



M, 



pi 



10 



-5 



0.01 



1/2 



M 



pi , 



H 



1/3 



r \i/2 



0.01 



; 



M, 



pi 



(3.36) 



(3.37) 



Hence, seeing tensors (r > 0.01) would put A close to the Planck scale, for couplings to hidden 
sectors that have order one non-Gaussianity (|/i| ~ H). 

Gaussian hidden sector. — If the hidden sector is Gaussian (|/Li| <C H), the constraint on high-scale 
physics comes from the non-linearity of the mixing operator in the regime of strong mixing. We use 
(3.26) to replace i>o in (2.12) and write (3.27) as 



_ Mm) ( H \ 

^NL- ^4/3 \^2^Ay' 



4/3 



(3.38) 
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where fi{M) = fi{M) d'^'^[M). The hmit M <^ H is particularly interesting, since the signal then 
only depends on the ratio of A and H. Solving for A, we obtain 

A = 5 X loM ^^^ J H . (3.39) 

Since Planck constrains equilateral non-Gaussianity at the level |/nl"^ | < 75, we conclude that 
A > 110 H. This limit does not depend on assumptions about the strength of the self-interactions in 
the hidden sector — in fact, S could have purely Gaussian correlations — but it does assume that the 
hidden sector field is light enough to contribute to the curvature perturbations: namely, we must 
have M <C p.^^ When there are nontrivial self-interactions in the hidden sector, the bound becomes 
stronger: see fig. 9. 

4 Supersymmetric Hidden Sectors 

We obtained the strongest constraints on high-scale physics, A > W^H, if the hidden sector field 

5 was light enough to be quantum- mechanically active during inflation (m < ^H) and had a 
cubic coupling obeying fj, > H. As we now explain, these conditions are very naturally met if 
supersymmetry is spontaneously broken during inflation [15]. On the other hand, we will also point 
out that if the hidden sector containing S is sequestered from the source of supersymmetry breaking, 
S can naturally have nearly Gaussian correlation functions. In this case, only our weaker — but 
universal — constraint, A > lO^H, applies. 

4.1 Local Non-Gaussianity from Generic Soft SUSY Breaking 

A supersymmetrized version of the inflationary theories discussed in this paper can be modeled 
in terms of three sets of superfields: a spurion X, an inflaton superfield $, and a hidden sector 
superfield S. The spurion X is taken to be the sole source of supersymmetry breaking. We assume 
that this SUSY breaking is communicated to $ and S only through gravity, i.e. via Planck-suppressed 
interactions. On the other hand, ^ and S can have direct couplings suppressed by a lower scale, 
such as the mixing term (2.5). The spontaneous breaking of SUSY by X induces soft terms in 
both the inflaton sector and the hidden sector. This situation is analogous to the generation of soft 
terms in the supersymmetric Standard Model through the spontaneous breaking of supergravity.^^ 
A complete characterization of SUSY-breaking soft terms in effective supergravity can be found 
in [33]. Here, we only review the basic conclusions relevant for our analysis. 

Generic soft terms. — We assume that SUSY is exact at high energies, but becomes spontaneously 
broken at low energies. In the low-energy EFT, the SUSY breaking is characterized by soft terms 
and the details of the breaking mechanism decouple. In the model of fig. 10, these soft terms can be 
computed in terms of A/" = 1 supergravity input [33] . For generality, we now allow both X and S 
to be collections of superfields, i.e. X* and S^, with z = 1, . . . , nx and I = 1, . . . , nj^. We take the 



^^Recalling that M = m^ /p and that the mixing is strong when p ^ H, we see that E can be light enough to 
contribute through strong mixing (m^ ^ p^) even ii m ^ H. In other words, fluctuations of E decouple from the 
curvature perturbations for m > p, rather than for m > |H as in the weak mixing case. 

^^In the application to the MSSM, E would stand for the 'visible' sector, while X is called the 'hidden' sector. Here, 
we have two hidden sectors, one that breaks SUSY (X) and one that receives non-Gaussian soft terms (E). 
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Figure 10: Couplings between the superfields in the supersymmetric effective theory of inflation. 

superpotential VF to be a general holomorphic function of X and S, and expand around S = 0: 

w{x, s) = w{x) + ^/i,j(x)s^s^ + ^y7ji(x)s^s-^s^ + • • • . (4.i) 

Famously, W is not renormalized at any order in perturbation theory [34], so that symmetry struc- 
tures imposed on the tree-level superpotential receive at most nonperturbatively small corrections. 
The Kahler potential K takes the following form 

K{X, X, S, S) = k{X, X) + Zjj{X, X)S^S^" + (^HijiX, X)S^S^ + h.A +■■■ . (4.2) 

All couplings in K receive perturbative corrections. The scalar potential for the spurion X is 

V{X, X) = kijF'P - ?,e^\W\^ , (4.3) 

where Kij = didjK and F^ = e^/'^K^^{diW + WdiK). We assume that the potential has a stable 
minimum with at least one non-zero F-term {F^) ^ 0, so that SUSY is spontaneously broken. The 
order parameter that measures the SUSY breaking is the gravitino mass, 

m3/2 = e<^)/2|M^((X))|. (4.4) 

Integrating out the spurion, one arrives at an effective theory for the field S, 

W^eff(S) = ^mS'S-^ + ^y^jiS^S^S^ , (4.5) 

where /x/j = e^^^^ifiu) + my^iHij) - {F'd,Hij) and Yjjl = e^^^l^{YijL). This induces the 
following soft SUSY breaking terms for S,^^ 

Vsofm = m]jT.'t^~ + Q^jjiS^S^S^ + ^S/jS^S-^ + h.c\ , (4.6) 

where 

"^/j = ^l/2^iJ ' "^IJL = FWiYijL , and Bjj = F^D^iij - 171^/2 IJ-IJ , (4.7) 



^For simplicity, we only show the result for Zjj = Sjj. The complete answer can be found in [33]. 
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with DiYjjL = diYjjL + \{diK)YjjL and Di^jj = di^jj -{- \{diK)iiij. We see that the generic scale 
of all soft terms is set by the gravitino mass: 



^]j ~ "^3/2 ' ^IJL ~ "ig/aY/jL , and Bij ~ m3/2Hij ■ (4.8) 

The final step is to relate 7713/2 to the scale of inflation, in order to compare the soft terms to 
the Hubble scale H. The precise relation is model-dependent, but reasonably general estimates are 
possible. We assume that instabilities of the moduli — for example, the decompactification instability 
of a string compactification — are controlled by barriers in the F-term potential. Comparing (4.3) 
and (4.4), we see that the size of the barriers is set by 7773,2 (i^i units with Mpi = 1). With the 
further assumption that the curvature of the moduli potential around its minimum is dictated by 
the scales appearing in W and K (rather than being much larger than these scales as a result of 
fine-tuned cancellations, cf. [35]), then the F-term energy can be at most of order the barrier height, 
KijF^F^ < few X 777? ,2- otherwise the moduli become destabilized. Since the F-term energy also sets 
the scale of inflation, we conclude that in generic configurations,^^ the scale of soft terms in SUSY 
inflation is the Hubble scale H. Thus, the masses and the cubic couplings of the hidden sector fields 
are quite generally at least of order H. Of course, mass terms and cubic couplings are allowed by 
supersymmetry, and the supersymmetric masses could be large enough to prevent S from fluctuating 
during inflation. When the supersymmetric contributions to the masses and couplings are small — for 
example, if E is a modulus, with a vanishing potential in the supersymmetric limit — then the soft 
terms control the dynamics. 

Comparing to our analysis of the non-Gaussianity in terms of the parameters 777 and //, we 
conclude that if the self-interactions of S are produced by generic soft supersymmetry breaking, 
then m ^ fi r^ H, and the dominant source of non-Gaussianity is self-interactions of S. The 
constraint on the scale of the dimension-five mixing operator is then A > W^H. 

4.2 Equilateral Non-Gaussianity from Sequestered SUSY Breaking 

The results of the previous section hold for generic soft symmetry breaking, for which all soft terms 
are of order the gravitino mass 7773/2. If instead the soft terms in the S sector are small compared 
to 7773/2, we say that S is sequestered from the source of supersymmetry breaking. Assuming further 
that S does not have large supersymmetric interactions — as above, S could be a modulus — we have 
jjL -^ H, and E has nearly Gaussian correlation functions. As we have shown, such Gaussian hidden 
sectors imprint equilateral or orthogonal non-Gaussianity in the curvature perturbations as a result 
of non-linearities in the mixing with #. In short, we conclude that the dominant channel for a sector 
of sequestered moduli to generate non-Gaussianity is by mixing with the inflaton. 

It is therefore worth asking whether sequestering of a hidden sector is plausible: is it guaranteed 
by reasonable symmetry assumptions, or in the context of an ultraviolet theory? Let us address 
these points in turn, in supergravity and in string theory. 

Sequestering in supergravity and in string theory. — Randall and Sundrum [38] originally observed 
that sequestering, in the form of vanishing tree-level soft terms, follows from the separability condi- 

^^This relationship was first emphasized in [.35]. Constructing natural models with 7713/2 ^ H remains a model- 
building challenge despite some dedicated efforts, cf. e.g. [:36, 37], which suggests that 7713/2 ~ H could be a relatively 
robust property of inflationary models in supergravity. 
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tions^" 

/(X,S) = /(X,X) + /(S,S), (4.9) 

W{X,^) = W{X) + W{^) , (4.10) 

where / = e^ '^. To understand whether these conditions are plausible, it is useful to recognize 
that sequestering amounts to the suppression of gravity-mediated interactions between the visible 
sector and the SUSY-breaking sector, beyond the degree of suppression determined by the Planck 
scale. In particular, the Kahler potential coupling 

d'^ex'^xT.^j: (4.11) 



Mp^i 



must have c <C 1. Because the coupling (4.11) is suppressed by the Planck scale, a proper justification 
for c <S 1 requires an ultraviolet completion of gravity. Unsurprisingly, the condition c ^ 1 is rather 
difficult to justify through symmetry arguments in four-dimensional effective quantum field theory.^^ 
For this reason we turn to extra-dimensional arguments for sequestering, in supergravity and in string 
theory. 

Before proceeding, we remark that the condition (4.10) might appear easy to justify, because if 
the tree-level superpotential is separable, then perturbative nonrenormalization [34] extends super- 
potential separability to all orders in perturbation theory. However, separability at the nonpertur- 
bative level does not generally follow, and so superpotential separability is a reasonable assumption 
only when nonperturbative superpotential terms can be neglected [40]. In fact, there are explicit 
examples^'^ in string theory in which nonperturbative separability is violated [42]. 

The original proposal for sequestering of supersymmetry breaking [38] relied on locality in an 
extra dimension: the conditions (4.9) and (4.10) were shown to hold when the supersymmetry- 
breaking sector and the 'visible''^^ sector communicated only via higher-dimensional gravity, with 
no other modes propagating in the bulk. Because a full ultraviolet completion is desired, one should 
determine whether the assumption of a barren extra dimension is justifiable in a compactification of 
string theory, i.e. does spatial separation imply sequestering in string compactifications? In general 
compactifications, spatial isolation alone does not suffice for sequestering [43], because light fields 
propagating across the compactification mediate interactions that induce large soft terms (cf. [44]). 
However, sequestering can be achieved if the visible sector and the hidden sector are separated along 
a warped extra dimension, e.g. if these sectors reside on D-branes at opposite ends of a warped 
throat region [45]: this is the gravity dual of conformal sequestering [41]. 

To summarize the above, the plausibility of sequestering should be determined in an ultraviolet- 
complete theory. In string theory, sequestered sectors can arise on D-branes in warped regions, 
though one must carefully check that nonperturbative superpotential couplings do not spoil seques- 
tering. 



^"Weaker conditions suffice to get partial suppression of the soft terms [39, 40]. 

^^The only complete four-dimensional argument [41], known as conformal sequestering, assumes that the hidden 
sector is strongly coupled and nearly superconformal over a large range of energies, and as such is dual to an extra 
dimension. 

^^ Concretely, a visible sector residing on D-branes can receive large soft terms through threshold corrections to a 
gaugino condensate superpotential generated on a distant stack of D7-branes [40, 42]. 

^^In the literature on sequestering, this is the visible sector containing the Standard Model, but for the present 
discussion it is actually the hidden sector containing S. 
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Properties of a sequestered hidden sector. — A natural question now is whether the particular struc- 
tures that ensure sequestering, particularly warping, can have any consequences beyond the smallness 
of the scalar soft terms. That is, do the sequestered scenarios arising in string constructions suggest 
different signatures from those described above in the $-S model? A general observation is that if 
warping is responsible for sequestering, then the dual approximately conformal field theory will be 
strongly coupled, with large anomalous dimensions. Non-Gaussianity from such a conformal hidden 
sector was characterized in [12], from the perspective that a conformal hidden sector is a reasonable 
possibility a priori. Here, we are pointing out in addition that in the context of spontaneously broken 
supersymmetry, one has two broad choices — generic soft breaking, leading to non-Gaussianity from 
self-interactions in the hidden sector, and sequestered breaking, with non-Gaussianity from mixing 
with the inflaton — and the latter scenario plausibly corresponds to a nearly conformal hidden sector. 
This provides further motivation for considering the theories described in [12]. 

5 Conclusions 

We have shown that the Planck limits on non-Gaussianity [3] imply constraints on hidden sector 
fields E coupled to the inflaton $ by nonrenormalizable operators suppressed by a scale A. We 
constructed a general EFT for ^ and E, assuming that ^ was invariant under a shift symmetry, 
but imposing no such restriction on S. The leading mixing between the two sectors comes from the 
dimension-five operator 

Non-Gaussianity in the curvature perturbations arises from two distinct sources: self-interactions 
in the hidden sector and nonlinear couplings between the two sectors. When the correlations of 
S have order-one non-Gaussianity — which is the case for natural-size cubic couplings — then hidden 
sector self-interactions dominate, and the curvature perturbations acquire non-Gaussianity of the 
local or equilateral type, depending on the mass of S. If instead S has a small cubic coupling and 
hence Gaussian correlations, the dominant non-Gaussianity in the curvature perturbations arises 
from nonlinearities in the S-$ coupling, and is of the equilateral or orthogonal type. Our simple 
two-field EFT therefore realizes all three bispectrum shapes probed by the Planck satellite. 

Although our analysis did not rely on assumptions about ultraviolet physics, it is worth pointing 
out that the effective theory described here is very natural from the perspective of string theory. 
The key ingredients are an inflaton field; one or more additional scalar fields that are light enough 
to fluctuate during inflation — for example, moduli that acquire masses m ^ H after supersymme- 
try breaking — and field-dependent kinetic couplings of the form (5.1). All three ingredients are 
commonplace in flux compactifications [11, 46]. 

Our results demonstrate that the Planck limits on non-Gaussianity constrain high-scale physics. 
Operators that mix a hidden sector involving light scalars with the inflaton must be suppressed by 
at least A > W^H. When the cubic couplings in the hidden sector are of order H, the limit is 
significantly stronger, A > W^H. Furthermore, a detection of primordial tensors [47], r > 0.01, 
would place severe constraints on non-Gaussian hidden sectors, by pushing the lower bound on A 
almost to the Planck scale: 

I I \ 1/3 , .1/9 



^^"Hir (m) ""■ (^^^> 
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The bispectrum results of Planck would then imply constraints on Planck-suppressed couplings 
to hidden sectors. Such constraints would provide a powerful selection principle for ultraviolet 
completions of inflation, because in known constructions — such as string compactifications in the 
supergravity regime — the effective cutoff scales, set by the Kaluza-Klein and string scales, are well 
below the Planck scale. Indeed, the presence of a cutoff that is parametrically below Mp\ is a 
hallmark of theoretically controllable ultraviolet completions of gravity. We find it noteworthy that 
a detection of primordial tensors would sharply constrain the existence of light hidden sector scalars 
in any such ultraviolet completion. 
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H^:i = ^^^^T^ and Hl^ = ,a\ (A.4) 



A Details of the Numerical Analysis 

In this appendix, we describe in detail the numerical analysis that produced the results of Section 3. 
The starting point is the Hamiltonian of the coupled ip-a system 

H = \^ + ^^{d.^f + \a^ + ^^{d,af + \m^a^ - p^a + ^^+/^a^ • (A.l) 

V J N- J 

V V 

Ha Hint 

The quadratic part, Hq — pipa, implies the equations of motion 

ip" + 2nip' + k'^^p = -pa[a' + 3na] , (A.2) 

a" + 271(7' +{k'^ + m'^a'^) a = pa ip' , (A.3) 

where Ti = a' /a and primes are derivatives with respect to conformal time r. The solutions to these 

equations determine the mode functions in the interaction picture. The effects of the two cubic 

interactions, 

r(i) _ 1 (5^V _^ ^(2) _ _3 

2 A 

are treated perturbatively. Throughout this appendix, we approximate the spacetime as de Sitter 

with scale factor 

a{T) = - , (A.5) 

r 

in units where H = 1. 

A.l Weak Mixing Approximation 

We refer the reader to [23] for a detailed discussion of the weak mixing regime {p <C H). In this 
limit the analysis simplifies because we can treat the mixing as a perturbative interaction. The 
interaction picture mode functions satisfy uncoupled equations of motion whose solutions are known 
analytically, at least in the quasi-de Sitter approximation. The bispectrum associated with the cubic 
self-interaction a^ gets converted to a bispectrum of curvature perturbations C, via three insertions of 
the quadratic mixing interaction (cf. fig. 5). In the in-in formalism, both the cubic interaction vertex 
and the insertions are integrated over time. The integrals have to be performed numerically. Results 
can be found in [23, 48]. Our treatment of the strong and intermediate mixing regimes {p > H) is 
new and will be described in the remainder of this appendix. 

A.2 Strong Mixing Approximation 

Parts of the strong mixing regime were studied analytically in [8]. Here, we develop a numerical 
treatment that covers the entire range of parameters. 

A. 2.1 Mode Functions 

In §2.3.2 and §3.2.2, we described the dynamics of the model in the strong mixing regime. We 
showed that the system reduces to a single degree of freedom whose evolution we can determine 
numerically. 
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Equation of motion. — For p ^ H, we can drop ip" + 27iip' on the l.h.s. of (A. 2) and a" + 2T-Lcr' on 
the l.h.s. of (A. 3). The mode function ^k{T) then satisfies 



.?-?(!+ ^'''>" 



''Jf't+^{{kr)' + m')n = 0. (A.6) 



where (p{u) satisfies 



It is convenient to rescale the time coordinate, u = kr/^J^, and write 

pi/4 
y^k{r) = -^-Viu) , (A.7) 

^"-^(l + ^^]^' + {u' + M)^ = 0, M="^. (A.8) 

We should think of ip as the solution corresponding to k = p = 1. The solutions for general k and p 
then simply follow from eq. (A.7). The solutions to eq. (A.8) have to be found numerically, except 
in the special case M = 0, where the answer can be written in terms of Hankel functions [8], 

^{u) = A{-uf/' K%iW) + B{-nf/' ^':;,{W) , M = . (A.9) 

Initial conditions. — In order to define the initial conditions, it is convenient to write q = a^ip. (This 
removes the friction term.) Moreover, at early times, /cr — )• — oo, we can expand the equation of 
motion in powers of the small ratio rr? jikT^'^ . At lowest non-trivial order, we find 

" V ' 

This equation has the WKB solution 

hm c?,(r) = ^^e-/-'=«<i^ (A.ll) 

fcr^-oo V'^fc(T) 



where oJkij) ^ t ~ 2~^ "^^^ constant ^fc = ^- is chosen so that the mode functions satisfy 

the Wronskian normalization g^g^' — Qlq'k = ^^ [8]- Eq. (A.ll) then implies 



t2 
1 



^3/2 ,U\Unn{-u) 



lim <f{u) = -= i-ur' 62" 62*^""^-"^ . (A.12) 

lt-5— oo ^2 

Wick-rotated solutions. — In order to improve the convergence of the bispectrum integral in the UV, 
we consider Wick-rotated mode functions. The WKB solution (A.12) is suppressed at early times if 
we choose the Wick rotation u — )• e'^^u. We define Wick-rotated mode functions as 

<pUu) = ^{e'^u) . (A.13) 

Since the mode functions are not known analytically, the Wick rotation has to be implemented at the 
level of the equation of motion. In particular, along the Wick-rotated integration contour, eq. (A.8) 
becomes 

^'M - - f 1 + / ... ) ^'Ju) - {u^ - iM) ^^{u) = , (A.14) 

u \ u^ — iM I ' 



with initial condition 



1 



lim ^^{u) = -^ {-uf' e-5^"^-*^>'^'-i-«»e-t*^ . (A.15) 

n-5--oo yJ2 
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A. 2. 2 Power Spectrum 

The superhorizon limit of the power spectrum of curvature perturbations is 






d^{M) 



H^ ( p\V 



^l\H 



where in the final equality we have restored explicit factors of H and defined 

d2(M) = 2|(^(0)|2 . 
A. 2. 3 Bispectra 

The cubic interactions in eq. (A. 4) induce a bispectrum for the inflaton fluctuations, 

{^kiVk2^k3) = C^-n-f B^{ki, k2, h) 5{ki + ^2 + fcs) • 
In what follows, we will use the notation 



'^w,k{v) 



1 



A;3/2 



• 'fwikv) . 



\ ^ iT- I n r\ Tvi r r^-y i~i /-i-r n /~vn I— I ' ' _ ^ ' ' 



Bispectrum from (dip) a. — The interaction H^^^ = ^ "' "^ gives 

A 



B^;\ki,k2,k3) = ^I''\ki,k2,k3) , 



(A.16) 

(A.17) 

(A.18) 
(A.19) 

(A.20) 



with 

I^'\ki,k2,ks) 



-Re 



ie 4 



dz; kl - kl - kl 



, V M - i{k^vY 

[VwM (0)<^^,fci {v)] Vpw,k2 {^)'P*wM (^)] \f^M i^WwM (^)] + 2 perms. 

Here, we have ignored the contribution from {ip)'^a, which is suppressed by powers of H/p. 
Bispectrum from, a^ . — The interaction H-^^^ = pa^ gives 



(A.21) 



B^^^iki,k2,k3)=pl^'\ki,k2,k3) 



(A.22) 



with 

X<''(A;i,A:2,/c3) = 12 Re 



ie^'i 



i. f^ dw V'li>,fcl(0)93^^fc^(f) 99^,^2(0)93;;; ;,2(v)99^,fc3(0)¥'^,fc3(^) 



V M-i{kivY M-i{k2vy M - iiksv^ 



B(^{ki,k2,k3) 






Full bispectrum. — The full bispectrum for the primordial curvature perturbations is 

^X3 - 

Using eq. (3.11), we compute the amplitude of the bispectrum 
Plots of the functions fi{M) and /2(M) are shown in fig. 6. 



where 



(A.23) 

(A.24) 

(A.25) 
(A.26) 
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A. 3 General Treatment 

Finally, we present a general analysis that is valid for all values of p, including the intermediate 
mixing regime p ^ H. Unlike the weak mixing regime, the interaction picture mode functions are 
not decoupled. Unlike the strong mixing regime, we cannot ignore half of the solutions. In this case, 
a numerical analysis seems unavoidable. The methodology that we will employ is a generalization 
of the treatment of [49] (see also [30]). 

A. 3.1 Quantization 

To discuss the quantization of the two-field system it is convenient to introduce the canonically 
normalized fields q^p = aip and q„ = aa. We combine these fields and their conjugate momenta into 
two-component vectors 

'"^ ^ --I^=f'rl- (A.27) 



dq' 
The fields mix via the following equations of motion 

€+(^^-^)l^ = +^U'^-l^-) ' (A-28) 

q'l + (fc2 + !I^) q^ = -^(^q'^ + i^^) . (A.29) 

In order to take into account the statistical independence of the two solutions we have to solve the 
equations twice. In the first run, the infiaton fiuctuation q^ starts in the Bunch-Davies (BD) vacuum, 
while the second field q^ is set to zero.^"^ We denote the corresponding solutions q^^' = {q\p , q^ Y' . 
In the second run, q^ starts in Bunch-Davies and q^ is set to zero. The corresponding mode functions 
are q^"' . We can combine the solutions from the two runs into a matrix with components Quj) = qj , 
with I, J G {(f, a}, i.e. 



Q={'tf)] ■ (A.30) 

\ qa' qa J 

The field operators are then written as^^ 

Qki'r) = Q{k,T)ak + h.c. , where afc = (afc ,0^ ) ■ (A.31) 

Imposing the following normalization conditions [49], 

(A.32) 
(A.33) 
(A.34) 



QQ*^- 


Q*Q^ = 


nn*^-n*n^ = o 


gn*^ - Q*n^ = a 


K\(^L^V] 


= 6^-^Sik - 



we get 

k') , (A.35) 

so that dj^ and (di , )^ are a set of independent annihilation and creation operators, respectively. 



^* Actually, this initial condition is only consistent in the limit p — >■ 0. For finite p, the equations of motion dictate 
a correction to the initial condition of qa- We treat this carefully below. 
^^In component form, this equation reads qi — Qi(j)a^''^ + (h.c). 
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A. 3. 2 Initial Conditions 

To determine the initial conditions for the coupled two-field system, we consider a WKB solution to 
the early time limit of the equations of motion, i.e. when u ^ p. 

WKB solutions. — At early times, we expect a; — )• /c + 0{t~^) and hence we have uj' — )• 0(t~^). 
We wish to solve the equations up to terms that are suppressed by r~^ (i.e. we will keep terms of 
order r~^). The WKB ansatz is 

Vw(r) 
q^ = ^^^e-'r^^'-^^'K (A.37) 

Plugging these into the equations of motion (A. 28) and (A. 29), and dropping terms of order r^^ or 
higher, we find 



{-u;'^ + k'^)q^fi + i-ujq„^o = , (A.38) 

P 



T 

{-Lo'^ + k^)q^fl-i(^uq^fi = . (A.39) 

Combining these two equations, we find 

uj^ -k^ = ±^oo . (A.40) 

r 

For each sign, this has two solutions corresponding to the positive and negative frequency solutions 
at early times. We focus on the positive frequency solutions 

u;±=k±f . (A.41) 

Substituting (A.41) into (A.38), we get a relation between q^^ and q^fl, 

where {(iitQ, Q^c^q) are the amplitudes corresponding to oj = uj±. We see that p has cancelled out, so 
that q^pfi and q^^ are of the same order. The WKB solutions therefore are 

v4/c 



V4/c 



where the overall normalization has been fixed by eq. (A. 34). To make contact with the discussion 
above, we define alternative basis functions 

q^^) = ^ ^^ . (A.46) 
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Using (A.43) and (A.44), we find 



if) {o-) 



1 ^_,,,/cos(nn(-A;r)) - sin(nn(-A;r)) 



p^o 1 _,w / 1 



Ilk Vsin(f ln(-A;r)) cos(^ ln(-A:r)) y V2/c 

(A.47) 
Tliese are tlie initial conditions we advertised above. Tliey are equivalent to the initial conditions 
defined in the (it) basis. In particular, the Bunch-Davies vacuum — defined as the state annihilated 
by both ayP' and aS"' — is also annihilated by a^^K In the following, we will find it slightly more 
convenient to work in the (it) basis. 

A. 3. 3 Mode Functions 

This time we have to solve the exact eqs. (A. 2) and (A. 3). Let us write 



V'fcW 



1 



(,5(u) and o"a;(t) 



1 



P/2 ^^^' "'^^^ /t3/2 

where u = kr. The equations of motion in a de Sitter background are 

3 



• o-(u) , 



-// 2 ,, p 

if — ^ +y^ = H — 
u u 



a' --a 



u 



~a" - 


2 , 


TV? 


--a' + 


1 + ^ 
U^ 




u 



a 



P ~i 
u 



(A.48) 

(A.49) 
(A.50) 



Initial conditions. — The two numerical runs are defined by the following set of initial conditions 

(A.51) 



^ -j(«+f ln(-«)) 



(+) : lim (^(+' =i-e-*^"+t 



(-) : 



It— > — oo 

lim if! 

u—^—oo 



lim c^(+) = ^e~*(''+2i''(~")) 
2 



w— >— oo 



(-) 



_i!fe-^("-5i'i(-")) 



f(-) 



u 



lim a^ ' = —e 

u—^—oo 2 



-J(M-f In(-u)) 



(A.52) 



Wick-rotated solutions. — This time we perform the Wick rotation u — )• iu. We define Wick-rotated 
mode functions as 

(pw{u) = (p{iu) and o"^(n) = a{iu). (A. 53) 

Eqs. (A.49) and (A.50) become 



^l 


-l^^- 


^w 


= +^ 

u 


., 3, 

U 


) 


(A.54) 


U 


V 21 

ra 


O'w 


p ~l 


(A.55) 




lim aW = ^e"-*^''^(-«)ef^ 

«— >•— oo 2 


(A.56) 


= -i-e-+^5 


in(-«)e-^ 


^^ 


lim 




«+ifln(-n)g-fp _ 


(A.57) 



with initial conditions 

(+) : lim ^'^' 

M— > — OO 

(-) : lim c^^-' 

M— >■— oo 2 n— >— oo z, 

A. 3. 4 Po^ver Spectrum 

The power spectrum of curvature perturbations is an incoherent superposition of the results of the 
two runs. 



2a2 



(27r)^A, 



C 



2 r 



|^'+'(0)|' + |^<-)(0)f 



(A.58) 
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A. 3. 5 Bispectra 

We define the following momentum-dependent mode functions: 

Wick contractions between the fields (p and a are 

r 



», 



'Pw,k{'>^l)'Pw,k{'"2) = {0\'Pw,k{^l)'Pw,k{'"'i)\^) = ^ 'Pw!k('"i)['Pw,k('"'2)] 



a=± 



.("), 



'Pw,k{vi)^w,k{'"2) = {0\f>w,k{Vi)^nj,k(.^2)\0) = Yj ¥'£fc(^l)[Vfc(^2)] 



(A.59) 

(A.60) 
(A.61) 



a=± 



Bispectrum from (d(p)^a. — The interaction H- 



(1) ^ 1 (8y)V 
int 2 A 



gives the bispectrum 



B^^\kuk2,h) = jI^'Kki,k2,k3) , 



(A.62) 



where X*^' is the sum of two terms, 1^2^ and I^g-^^^^, which arise from the interactions generated by 
the ip'^a and {dnp)'^a terms, respectively. These two terms read 



1^2^{ki,k2,k3) = 2Re 



dv,[ 



['fwMi^)'P'wM('")] 



X [0w,k2iO)'f'w,k2('")]['Pw,k:iiO)^y^Miv)] +2 perms. 



(A.63) 



2^0,^)2.(^1,^2,^3) 



-Re 



^y 

-2 \h - h- ^i) 



X [0wMi^)'PwMi'^)] [0w,k2i'^)^w,k2i'^)] [^■w,k3i'^)^wMi'>^)] + 2 perms. 

(A.64) 

Bispectrum from, a^ . — The interaction H^^l = iia^ gives 

B^^\ki,k2M) = k^^''\kiM,k^) , (A.65) 

with 

/o df ' ' ' ' ' ' 

-oo ^ 

(A.66) 



Full bispectrum,. — The full bispectrum of the primordial curvature perturbations is 



B^{ki,k2,k3 



H 

6n 



^(1) , ^(2) 



Using eq. (3.11), we get 



/nl = fi{m,p) ■ — + f2{m,p) 



1 fi 



27rA^ H 



(A.67) 



(A.68) 
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0.01 0.1 1 10 100 

Figure 11: Numorical computation of fi{m, p) (dashed) and /2(m,p) (solid) for m = H. 



where 



fi{m,p) = 



5 XW(fc,fc,fc) 
18 Pj(fc) 



and 



f2{m,p) 



-27rA„ 



5 I'-^\k,k,k) 
18 



P'ik) 



(A.69) 



IR divergences. — The factorized form of the in-in formula is not well-suited to deal with the IR 
behavior of the integral, i.e. the behavior as f — )• 0. This leads to a spurious IR divergence when the 
integral is computed numerically. In order to make the IR convergence manifest, one needs to use 
the so-called commutator form, which schematically reads 



{^' 



dT{[Hint,^^] 



(A.70) 



However, the Wick rotation is no longer valid in the commutator form and therefore one cannot 
achieve convergence in the UV. Hence, one way to make the integral converge both in the UV and 
the IR is to split it in two parts [23] 



(^' 



5,3 \ 



2Re 



-oo(l— je) 



+ / dr([i?int(r),^3]^. 



(A.71) 



where the UV part (r < Tc) is computed in the factorized form and then Wick rotated, while the 
IR part (r > Tc) is computed in the commutator form. This solves the problem of spurious IR 
divergences. Naturally, the final result is independent of the choice of cutoff Tc used to distinguish 
the UV and IR parts of the integral. 
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B Details of the Effective Theory 

The operator (2.5) that mixes the inflaton sector ^ and the hidden sector S contains a tadpole for S, 

/:mixD(/9^0)S . (B.l) 

In Section 2, we assumed that a potential V{Ti) stabilizes the field at Eq and then studied the 
phenomenology of fluctuations cr = S — Sq. In this appendix, we will discuss this problem in more 
detail. Because ^o ^ H^i we expect the minimum Eq to be displaced significantly from the origin. 
The effective couplings of fluctuations around Eq can therefore be quite different from those near the 
origin. In the cases of greatest phenomenological interest, all dimensionful couplings in the effective 
action for a were of order the Hubble scale H. Here, we want to explore whether this structure is 
natural after including the vev for E. We will describe two different vantage points: the effective 
theory of the background [18] and the effective theory of fluctuations [7]. For concreteness, we will 
limit the discussion to the weak mixing regime. 

B.l Effective Theory of the Background 

We start with a simple model in which a single scale controls both the mixing term and the hidden 
sector self- interactions. For weak mixing, this model will fail to produce measurable non-Gaussianity, 
but the way in which it fails will be instructive. In particular, it will motivate models in which the 
mixing interaction and the hidden sector self-interactions are controlled by two distinct scales. 

One-scale models. — Consider the canonical slow-roll Lagrangian 

U = -^id^f-V{^). (B.2) 

We include a tree- level coupling to a free field E, 

1,^^,2 l(a$)2 



Cs + C 



mix 



2(5S)^-2^S- (B-3) 



We have chosen the classical potential for E to vanish, so that E is protected by a shift symmetry 
that is broken only by the higher-dimensional mixing term. All corrections to the potential for 
E therefore vanish as A — )• oo. Let us assume that this effective description holds up to a cutoff 
scale Aj,. The one-loop effective action for E is then 



Cj:+Cr 



i-ioop" 2^"^^' 2 A ^ 167r2A2^ 



2 



l + c^l+c,^ + 



(B.4) 



We have cancelled the one-loop tadpole, since we are assuming that (E) = in the vacuum. In order 
to induced a mass of at most m?, we require that A^ < IGvr^m^A^. Using this mass term, tti^E^, to 
stabilize the tadpole for E, we find 

For strong mixing, we require /) S> M = m'^ /p, which implies that Eq ^ A. As a result, in the case 
of strong mixing, the EFT description breaks down between E = and E = Eq . We will discuss the 
EFT around Eq in the next subsection. 
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The case of weak mixing {p < H) is much less constrained: we are free to take p < m < H, so 
that Sq < A. When this holds, higher-order terms in the potential are suppressed by Sq/A and are 
therefore not important. The structure of this effective theory is appealing: the small mass for E is 
explained by the approximate shift symmetry for S (assuming an appropriate UV completion above 
Aj,). Furthermore, we have found that the vev, Sq, is under control: higher-dimension operators are 
suppressed by powers of Sq/A. 

Unfortunately, for weak mixing this setup is too restrictive to allow for measurable non-Gaussianity. 
To see this we note that the cubic coupling in (B.4) is // ~ H'^/A. Eq. (3.17) then implies 

where we have used (p/H)'^ < c^^{m) as required for perturbative control of the power spec- 
trum (3.16). We see that the amplitude is parametrically suppressed and unobservable in present 
and future CMB experiments. Including a tree-level potential V(T,) = //S^, with p, ^ H, does not 
improve the situation. At Sq, this would induce a mass term m'^ ~ H'Eq. To keep the field light, we 
have to require T,o ^ H (or accept fine-tuning against the bare mass.). But eq. (B.5) then implies 
p < {2ttA^)H and we get 

U < i27rA^ff{m) ~ 10"^ f{m) < 1 . (B.7) 

We conclude that the one-scale model, although natural, is too weakly coupled to produce a signifi- 
cant three-point function. The same applies to the four-point function, 

— <")W'(f)^(l)'-(»)(l)'<^^^. (-' 

where we have used the quartic coupling in eq. (B.4). Although we have addressed the small size of 
the mass term and cancelled the tadpole, we have also suppressed the interactions that generate the 
non-Gaussian correlations. 

Two-scale models. — For weak mixing, in order to produce measurable non-Gaussianity, we need an 
additional cubic interaction beyond the one induced by the mixing term. The basic problem we just 
ran into can be traced back to the fact that we associated the same scale to the breaking of the 
shift symmetry and to the strength of the cubic (or quartic) interaction. The obvious strategy to 
produce a measurable signal is to separate these two scales. The examples that we will present are 
meant as simple existence proofs of natural theories with large non-Gaussianities, rather than as a 
comprehensive study of the range of possibilities. 

The most straightforward modification to the effective Lagrangian in (B.3) is to introduce ad- 
ditional interaction terms that do not violate the shift symmetry for E. These interaction terms 
can then be large without any risk of introducing large corrections to the mass of S. As a concrete 
example, we will consider 

After including the one-loop corrections (B.4), the tadpole is stabilized as before. At Sq, we have 
the quartic interaction 

Ant = Y^ida^ . (B.IO) 
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The associated trispectrum amplitude is 

''-=«<'">-(24p'(f)'(^)*' "*■"' 

For A ^ A, this can be significantly boosted relative to (B.8). An observable trispectrum can arise 
if A < (<i>o)^ • The scale A is not consistent with the effective theory for the background at a scale 
of order (^q)^' ^. Nevertheless, couplings of this type are allowed within the effective theory for the 
fiuctuations, provided that A> H = 0.02 (^q)^- Furthermore, the dynamics that lead to this large 
interaction are confined to the S sector and do not affect the evolution of $o(i)- 

The above model illustrates how to introduce additional interactions for S without affecting the 
mass of the field. In principle, the same logic can be used to produce a large bispectrum. A simple 
operator that achieves this is 

^^d,^d^^{d^f ^ i^^aidaf . (B.12) 

For small enough A, this allows measurable non-Gaussianity. The operator also contains additional 
mixing terms that transfer the fluctuations to the visible sector. It would be interesting to explore 
these model-building considerations further, but this is beyond the scope of the present paper. 

B.2 Eflfective Theory of the Fluctuations 

In expanding V{T,) around So, we have assumed a single EFT for the background that is valid 

around both S = and S = Sq. Imposing such a constraint is unnecessarily strong, as the a field 

only probes a small region around Eg. By contrast, in the effective theory of inflation [7], one writes 

down the theory for the fluctuations directly. The only remnants of the background solution are the 

(time-dependent) couplings in the action. Specifically, the Lagrangian for the fluctuations is given 

by 

1 .1 .If. 1 .1 

a-V{a). (B.13) 



£ = _i,a,)^_i(8,)^_i 



^o^--{d<ff 



By writing the mixing term as [^Qip — 2(c?'/')^]o", we have explicitly^^ removed the tadpole for a, 
which is stabilized at a = for all time. As in the main text, we take V{a) = i^m?a'^ + /uo"^ + • • • . 
It is the parameters A, m'^ and /x that are constrained by Planck's limits on non-Gaussianity. 

In the main text, we have shown that the mass and the cubic coupling receive the one-loop 
corrections 

5m2 ~ ^ and Sfi ^ ^ . (B.14) 

A2 A3 

In the case of weak mixing, m^ < H^ and p , fi < H, these are small corrections when A^ ~ ($0)^' ^. 
For strong mixing, p ^ H, these corrections are more dangerous but still produce a natural EFT 
when A* ~ p < (<i>o)^ • One may be tempted to conclude that there is no flne tuning in the EFT 
of the fluctuations. On the other hand, this does not seem to match our experience in the case of 
slow-roll inflation with a polynomial potential for S. The resolution is that, for slow-roll inflation. 



^^This is more transparent in the language of the EFT of inflation [7], where this mixing term is written as /Imix ~ 
— $o[c'(i + 7r)^ + l]o" with tv = ifi/^. The tadpole had to be cancelled by hand by including the +1, otherwise the action 
would break diffeomorphism invariance explicitly. 
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the interactions in the effective Lagrangian are negUgible and we are free to take A* ~ A. Under 
such circumstances, we would find Sm? ~ A^ ^ $o- We conclude that without a new scale below 
A, writing a small mass for a is not natural (for weak or strong mixing). 

B.3 Summary 

In the previous two subsections, we described the naturalness of the EFT from the point of view of 
the background and the fluctuations. In terms of the background, we found that our EFT around 
E = can break down before we reach S = Sq. On the other hand, the theory of the fluctuations 
is well-defined and natural, provided that the theory is cut off at A^ < ($o)^ • 

To make contact between the two descriptions, let us consider an effective theory for the back- 
ground that is well-defined for a region of size A^ < A around Sq. For example, near Sq we might 
have 

F(S) = m4-"S" + m2(S - So)2 + Mo(S - ^of + • • • , (B.15) 

where 1 < a < 2. The first term in the potential is responsible for canceling the tadpole at S = Sq. 
Expanding the leading term around S = Sq, we find that the corrections to the mass and to the 
cubic coupling are given by 

Sm ~ — — and 6fi ^^ — k- . (B.16) 

Taking So ^ ^q/H ensures that these corrections are negligible. However, for the loop corrections 
to be small, our "effective theory" must be cut off at Aj, r^ {^)q ^ ^o- Therefore, this effective 
theory does not include S = 0. 

In conclusion, there exist effective theories for the fluctuations where the mass for a is naturally 
small, yet there is still a large non-Gaussian signal. However, a single weakly coupled effective theory 
does not include both the vacuum and inflationary values of S. 
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